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We explain the excess of the antiproton fraction recently reported by the AMS-02 
experiment by considering collisions between cosmic-ray protons accelerated by a local 
supernova remnant (SNR) and the surrounding dense cloud. The same u pp collisions” 
provide the right ratio of daughter particles to fit the observed positron excess simulta¬ 
neously in the natural model parameters. The supernova happened in relatively lower 
metalicity than the major cosmic-ray sources. The cutoff energy of electrons marks the 
supernova age of ~ 10 5 years, while the antiproton excess may extend to higher energy. 
Both antiproton and positron fluxes are completely consistent with our predictions in 
Ref. [4], 


1. Introduction 

Recently the Alpha Magnetic Spectrometer (AMS-02) on the International Space Station has 
reported that the antiproton to proton ratio stays constant from 20 GeV to 450 GeV kinetic 
energy [1], This behavior cannot be explained by the secondary antiprotons from collisions 
of ordinary cosmic rays with interstellar medium [e.g., 2], This suggests a new source such 
as astrophysical accelerators and annihilating or decaying dark matter, although there are 
still uncertainties in the background modeling [3]. 

The excess of antiprotons looks surprisingly similar to what we predicted [4] when the 
PAMELA experiment detected the positron excess [5] and the Fermi, HESS, and ATIC/PPB- 
BETS experiments observed the electron anomaly [6-8]. We considered recent supernova 
explosions in a dense gas cloud (DC) near the Earth. The antiprotons and positrons are 
produced as secondaries by the pp collisions between cosmic-ray protons accelerated by 
the supernova remnant (SNR) and target protons in a DC which surrounds the SNRs [4], 
Since the fundamental process determines the branching fraction, the positron excess should 
accompany the antiproton excess. 

There are several variants of such a hadronic model, e.g., the reacceleration of secondaries 
by SNR shocks [9] or the non-standard propagation that increases secondaries from ordinary 
cosmic-ray collisions with interstellar matter [10-12], Since the element ratio is the same as 
the ordinary cosmic rays in these models, the ratio of secondaries (e.g., Li, Be, B) to primaries 
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(C, N, O) must also rise with energy beyond ~ 100 GeV, which is not observed yet [13-15]. 
In contrast to these models, our model can accommodate the observations as shown below. 

Before AMS-02, the antiproton observations were consistent with the secondary back¬ 
ground [16]. Thus the leading models for the positron excess are leptonic, such as pulsars 
and leptophillic dark matter [e.g., 17-20]. The pulsars cannot usually explain the antiproton 
excess. On the other hand, the dark matter for the positron excess is now severely con¬ 
strained by other messengers such as gamma-rays and cosmic microwave background [e.g., 
21, 22], and hence we may need fine tunings in dark matter models to reproduce both the 
antiproton and positron excesses [3, 23-25]. 

Following Occam’s razor, we reexamine our nearby SNR model and simultaneously fit 
the antiproton fraction, positron fraction, and total electron and positron flux in light of 
new AMS-02 data. In particular the pp collisions give the correct branching fraction for the 
observed positron to antiproton ratio. Throughout this paper we adopt units of c = h = 
ks = 1 . 


2. Supernova explosions in a Dense Cloud 

Here we consider supernova explosions which occurred around ~ 10 5 — 10 6 years ago in a 
DC. We assume the DC is located at around ~ 100-200 pc away from the Earth like the 
progenitor DCs that produced the Local Bubble (LB) or Loop I. In general a massive star 
tends to be born in a giant DC [26] which explodes as a supernova. In this paper, we assume 
that the Giant DC is ionized and the temperature is approximately 10 4 K at the time of 
its explosion [27]. The shock of an SNR accelerates protons, which produce copious energetic 
mesons (pions and kaons, etc.) and baryons (antiproton, proton, antineutron, neutron, etc.) 
through the pp collisions in the surrounding DC. The mesons further decay into energetic 
positrons, electrons, gamma-rays, and neutrinos. In total those local secondary particles can 
be observed at the Earth as cosmic-rays in addition to the standard background components. 

The energy spectrum of the accelerated protons is parametrized by 
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where s is the spectral index. The age of SNR, i age , approximately determines the maximum 
energy [28], 


Emax ’ p ~ 2 X 1Q2 ^ 8 ^ 10 ^ lCLyr 
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where the shock velocity, v s , is v St $ = u s /10 8 cm s -1 ~ O(l)- Bd is the downstream magnetic 
field. We take the minimum energy of the protons to be its rest mass. We assume that the 
supernova explodes at the center of a DC for simplicity. In addition, we also assume that 
the acceleration stops when the Mach number of the shock decreases to 7 [4], We define this 
time as the acceleration time t acc = f age , and the energy spectrum at this time is given by 
s ~ 2 and E maXjP ~ 120 TeV [4], The SNR continues to expand even at i age > f acc . 

The radius is 50 pc at t age = 5 x 10 5 yr. Since it is comparable to the size of a giant DC, 
i?DC) [29] and the initial energy of the ejecta from the supernovae is larger than the binding 
energy of a DC, the cloud would be destroyed around this time. Until it is destroyed, the DC 
is illuminated by the accelerated protons from the inside with the spectrum of Eq. (1) given 
at t age ~ t acc - The duration of the exposure, t pp , could be approximated by the time elapsing 


2/8 







from the explosion of the supernovae to the destruction of the DC because the timescale t acc 
is shorter than 5 x 10 5 yr. 

After the destruction of the DC, the produced charged particles such as p, p, e + , or e~ 
propagate through diffusion processes and reach to the Earth. Since we assume that the 
DC has already been destroyed well before the present epoch, there are some differences in 
arrival times between those charged particles and massless neutral particles such as photons. 
It should be a reasonable assumption that we would not detect any photon and neutrino 
signals from the DC f 0 5-6 

years after the destruction of the DC. 

We have calculated spectra of those daughter particles through the pp collisions by per¬ 
forming the PYTHIA Monte-Carlo event generator [30] (See [28] for the details). Then we 
solve the diffusion equation of the charged particle “i” (i runs p, p, e + , and e~, ), 

Of- 0 

= K(ei)Afi + -q— [B(ei)fi\ + Q(ei) (3) 

where f t (t, x, £j) is the distribution function of an i particle, and £* = £//GeV with Ej being 
the energy of the i particle. The flux is given by 


$i(t,X,£i) =—fi. (4) 

We adopt a diffusion model 08-005 given in [31] with the diffusion coefficient, 

K (e e ) =K 0 (l+ 3GeV ) , (5) 

with A’o = 2 x 10 28 cm 2 s _1 and 5 = 0.42 [24, 32, 33]. The cooling rate through the syn¬ 
chrotron emission and the inverse Compton scattering is collectively parametrized to 
be [34] 


B{e e 


~ lO-^s- 1 ^ 




+ 0.9 


( 6 ) 


where B ^ft is the magnetic field outside the DC. This set of the parameters approximately 
corresponds to the MED model of the cosmic-ray propagation [35]. 

If we assume that the timescale of the production is shorter than that of the diffusion, 
~ d?/(Kc ) with d being the distance to the source, and the source of the daughter particles 
is spatially localized sufficiently, we can use the known analytical solution in [36]. When the 
shape of the source spectrum is a power-law with an index a to be 


Q = Q 0 s- a 5{x)8(t), 


(7) 


then the solution is given by 
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where £ cut = e^/Bidm, and the diffusion length is represented by 


/ 1 - (1 - -£«- 
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(9) 


d means the effective distance to the source by spatially averaging the distance to the volume 
element of the source, and we assume a ~ s. We approximately have 

d 2 rii 
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Fig. 1 Antiproton fraction fitted to the data. The data points are taken from [1] for AMS- 
02, and from [16] for PAMELA. The dotted line is plotted only by using the background 
flux [37]. The shadow region represents the uncertainties of the background flux among the 
propagation models shown in [1], Cosmic rays below an energy < lOGeV are affected by 
the solar modulation. We choose the background line and its uncertainty band only for a 
demonstration purpose. This choice is not essential for our conclusion (See the text about 
Fig. 3). 


with V s the source volume where 


d 2 rii 

dtdEi 



( 11 ) 


The differential cross section of the “/’-mode for the production of the i particle is rep¬ 
resented to be dcrj(Ep, Ei)/dEi with the multiplicity into the j-mode, gj = gj(E p ,Ei). 
v p = Vp(Ep) is the velocity of the primary proton. We also consider the free neutron 
(antineutron) decay for the electron (positron) production process, which is not included 
in the original version of PYTHIA. The initial proton spectrum can be obtained by a 
normalization to satisfy 

V s J dEp-j-j^r = Ftot.p. (12) 

For the local propagation of protons and antiprotons, their cooling is negligible unlike elec¬ 
trons and positrons. Additionally we can omit annihilations of antiprotons through scattering 
off the background protons because the scattering rate is small. We can also omit convection 
by interstellar turbulence within the galaxy. An analytical solution for the proton and the 
antiproton is also given by the same equation as Eq. (8) with a limit of s p /e cu t = 0. 
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Fig. 2 (a) Positron fraction (solid line), which includes the electrons and positrons coming 

from the DC and background electrons (dotted line, for example see Refs. [31, 34]). Filled 
circles correspond to the AMS-02 data [1, 38, 39] and PAMELA data [5] (b) Total electron 
and positron flux (solid line). The flux of the electrons and positrons created only in the 
DC (background) is plotted by the dashed (dotted) line. Observational data by AMS-02, 
Fermi, HESS, BETS, PPB-BETS, and ATIC2 [6-8, 40] are also plotted. The shadow region 
represents the uncertainty of the HESS data. 


3. Antiproton and positron fittings 

In Fig. 1, we plot the antiproton fraction at the Earth in our model (See also a similar model 
named “model B” given in Ref. [4]). For the background flux, we adopted the 20% smaller 
value of the mean value shown in [37]. Here, the radius of a spherical DC, 7?dc = 44) pc is 
adopted. The target proton density is set to be no = 50 cm -3 . The spectral index s = 2.15 
and the maximum energy, E max = 100 TeV, are assumed. We take the duration of the pp 
collision to be t pp = 2 x 10 5 yr. The total energy of the accelerated protons is assumed to be 
Ftot,p = 2.6 x 10 5 ° erg. The distance to the front of the DC is set to be d = 200 pc. About 
the diffusion time of e~ and e + , t^ifi = 2 x 10 5 yr is adopted. We take the magnetic field 
outside the DC to be B^m = 3 /iG (See [4] for the further details). 

In Fig. 2, we also plot the positron fraction and the total e _ +e + flux. It is remarkable that 
we can automatically fit the observational data of both the positron fraction and the total 
e~ + e + flux by using the same set of the parameters [4], Here the cooling cutoff energy is 
approximately given by e cut = el/Bt^m ~ 1 TeV (idiff/2 x 10 5 yrs)" 1 . 

The positron fraction rises at higher energies than that of the antiproton fraction (Fig. 3), 
because the spectral index of the background antiproton is harder than that of the back¬ 
ground positron. This comes from a difference between their cooling processes. Only for 
background positrons and electrons the cooling is effective in the current situation. 
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In Fig. 3, we plot the positron to antiproton ratio as a function of the rigidity. Here the 
local components represent the contribution of the nearby SNRs produced only by the pp 
collisions. From this figure, we find that both of the positron and the antiproton can be 
consistently fitted only by adding astrophysical local contributions produced from the same 
pp collision sources. 
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Fig. 3 Positron to antiproton ratio as a function of the rigidity with adding the local 
components produced by the pp collisions occurred at SNRs near the Earth. The thick solid 
line represents the case of the total flux. From the upper right to the lower left, we plot the 
flux ratios of 1) the one at the source (without cooling), 2) only the local components, 3) the 
total of the local and the background components, and 4) only the background components. 
The observational data reported by AMS-02 are also plotted. 


4. Conclusion 

We have discussed the anomaly of the antiproton fraction recently-reported by the AMS-02 
experiment. By considering the same origin of the pp collisions between cosmic-ray protons 
accelerated by SNRs and a dense cloud which surrounds the SNRs, we can fit the data of 
the observed antiproton and positron simultaneously in the natural model parameters. The 
observed fluxes of both antiprotons and positrons are consistent with our predictions shown 
in Ref. [4], 

Regardless of the model details, the ratio of antiproton to positron is essentially determined 
by the fundamental branching fraction into each mode of the pp collisions. Thus the observed 
antiproton excess should entail the positron excess, and vice versa. This does not depend on 
the propagation model since both antiparticles propagate in a similar way below the cooling 
cutoff energy TeY. 


6/8 








The cutoff energy of e cooling marks the supernova age of 10 5 years [19, 41], while 
we also expect a e + cutoff. The trans-TeV energy will be probed by the future CALET, 
DAMPE and CTA experiments [42, 45]. An anisotropy of the arrival direction is also a unique 
signature, e.g., [43]. We may estimate the amplitude of anisotropy as 5 e ~ 3d/2cidiff ~ 0.5%, 
which is below the upper limits by Fermi observations [44]. 

The boron to carbon ratio as well as the Li to carbon ratio have no clear excesses [1]. 
This suggests that the carbon fraction of the excess-making cosmic rays is smaller than that 
of the ordinary cosmic rays. In general the supernovae in the DC would not be the main 
channel of cosmic-ray production. Most of cosmic rays above ~ 30 GeV may be produced 
in chemically enriched regions, such as superbubbles, as implied by the hard spectrum of 
cosmic-ray helium [46]. Or the carbon abundance of the destroyed DC might happen to be 
lower than the Galactic average [4], 

We should be careful about the background systematics. In particular the propagation 
uncertainties yield the largest errors [3, 47]. However, in the energy region above ~ 100 GeV 
where the background contribution is small, the observed positron to antiproton ratio is very 
close to the branching fraction of the pp collisions (source components in Fig. 3). This fact 
is free from the background choice and partially supports our model. 
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Note added 

While finalizing this manuscript, Ref. [48] appeared which has some overlaps with this work. 
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